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p^i ■ Abstract. We consider the linear Dirac operator with a (— l)-homogeneous locally periodic 

5_l ■ potential that varies with respect to a small parameter. Using the notation of G-convergence 

for positive self-adjoint operators in Hilbert spaces we prove G-compactness in the strong re- 
solvent sense for families of projections of Dirac operators. We also prove convergence of the 
corresponding point spectrum in the spectral gap. 
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^3 ■ 1. Introduction 

S . 

In the present work we study the asymptotic behavior of Dirac operators J%h with respect to 

a parameter h € N as h — > oo. We consider Dirac operators J$?h = H + Vh on L 2 (M. 3 ; C ), where 

H = Ho + W + 1 is a shifted Dirac operator. The operators Ho, W, I, and Vh are respectively 

the free Dirac operator, the Coulomb potential, the 4x4 identity matrix, and a perturbation 

CN ! to H. We will study the asymptotic behavior of J% and of the eigenvalues in the gap of the 
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continuous spectrum with respect to the perturbation parameter h. 



G-convergence theory which deals with convergence of operators, is well-known for its ap- 
plications to homogenization of partial differential equations. The concept was introduced in 
the late 60's by De Giorgi and Spagnolo [5l [91 [TO] for linear elliptic and parabolic problems 
with symmetric coefficients matrices. Later on it was extended to the non-symmetric case by 
Murat and Tartar [3 [T2| [13] under the name of H-convergence. A detailed exposition of G- 
convergence for positive self-adjoint operators is found in Dal Maso [3]. In the present work we 
will base a lot of our framework on the results in Chapter 12 and 13 in |3j. The Dirac operator 
is unbounded both from above and below. This means that the theory of G-convergence for 
positive self-adjoint operators is not directly applicable to Dirac operators. In this work we 
study self-adjoint projections of Dirac operators which satisfy the positivity so that the theory 
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of G-convergence becomes applicable. 

We will consider periodic perturbations, i.e. we will assume that the potential V^ is a periodic 
function with respect to some regular lattice in R . We are then interested in the asymptotic 
behaviour of shifted perturbed Dirac operators J#h. This yields homogenization problems for 
the evolution equation associated with the Dirac operator J%?h 

{ Ufc(-,0)=u° 
and the corresponding eigenvalue problem 

jfe h u h (x) = \ h u h {x). 

The paper is arranged as follows: In Section 2 we provide the reader with basic preliminaries 
on Dirac operators, G-convergence and on the concepts needed from spectral theory. In Section 

3 we present and prove the main results. 

2. Preliminaries 

Let A be a linear operator on a Hilbert space. By R(A), D(^4), and N(A) we mean the range, 
domain, and null-space of A respectively. 

2.1. Dirac Operator. We recall some basic facts regarding the Dirac operator. For more 
details we refer to the monographs |14| . |15j and |17j . 

Let X and y denote the Hilbert spaces ^(M^C 4 ) and L 2 (R 3 ;C 4 ), respectively. The free 
Dirac evolution equation reads 

d 

(1) ih—u(t, x) = H u(i, x) , 

where Ho : ^ — > V is the free Dirac operator with domain ~D(Hq) = X, which acts on the 
four-component vector u. It is a first order linear hyperbolic partial differential equation. The 
free Dirac operator Ho has the form 

(2) H = -ihca ■ V + mc 2 (3 . 

Here a ■ V = X]j=i a *^"' ^ ^ s the Planck constant divided by 2-7T, the constant c is the speed of 
light, m is the particle rest mass and a = (a±, «2, 03) and j3 are the 4x4 Dirac matrices given 
by 

a-: = \ I and B = 1 . 

U ) I -I ) 
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Here I and are the 2x2 unity and zero matrices, respectively, and the <Tj's are the 2x2 Pauli 
matrices 

0-i=| , 02 = and 03 = 

V i ° J V i ° J V ° - 1 J 

Note that a separation of variables in (JTJ) yields the Dirac eigenvalue problem 

(3) Hqu(x) = Xu(x) , 

where u(x) is the spatial part of the wave function u(x, t) and A is the total energy of the 
particle. The free Dirac operator Hq is essentially self-adjoint on Cg°(M 3 ;C 4 ) and self-adjoint 
on X. Moreover, its spectrum, a (Ho), is purely absolutely continuous (i.e. its spectral measure 
is absolutely continuous with respect to the Lebesgue measure) and given by 

0-(Ho) = (— oo, — mc ] U [mc , +oo) . 

Ho describes the motion of an electron that moves freely without external force. Let us now 
introduce an external field given by a 4x4 matrix- valued function W, 

W(x) = W lJ (x) i,j = 1,2,3,4. 

It acts as a multiplication operator in L 2 (M 3 ;C 4 ), thus the free Dirac operator with additional 
field W is of the form 

(4) H = H + W . 

The operator H is essentially self-adjoint on Co°(M 3 ;C 4 ) and self-adjoint on the Sobolev space 



X provided that W is Hermitian and satisfies the following estimate (see e.g. [1 

(5) \Wij(x)\ <a^- + b, Vi£t 3 \{0} i,j = 1,2,3,4, 

2\x\ 

the constant c is the speed of light, a < 1, and b > 0. From now on we let W(x) be the 
Coulomb potential W(x) = —I, where Z is the elctric charge number (without ambiguity, I is 
usually dropped from the Coulomb term for simplicity). The spectrum of the Dirac operator 
with Coulomb potential is given by 

0(H) = (-co, -mc 2 ] U {A fc } fceN U [mc 2 , +oo), 

where {A fc }fcgN is a discrete sequence of eigenvalues in the "gap" and the remaining part of the 
spectrum is the continuous part o-(Ho). 

In the present paper we consider a parameter-dependent perturbation added to the Dirac 
operator with Coulomb potential. The purpose is to investigate the asymptotic behavior of the 
corresponding eigenvalues in the gap and the convergence properties. To this end we introduce 
a 4x4 matrix- valued function Vh = Vh(x) and define the operator J#h as 

(6) Jtfh = H + V h . 
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We recall that a function F is called homogeneous of degree p if for any nonzero scalar a, 
F(ax) = a p F(x). The next theorem is of profound importance for the present work. 

Theorem 1. Let W be Hermitian and satisfy the bound ^ above. Further, for any fixed h £ N, 
let Vh be a measurable (— 1) -homogeneous Hermitian 4x4 matrix-valued function with entries 
in Lf (M. 3 ), p > 3. Then J%h is essentially self-adjoint on C^°(R 3 ;C 4 ) and self-adjoint on X. 
Moreover 

o-(J^h) = (-oo, -mc 2 ] U {A^} fceN U [mc 2 , +oo), 

where {A^}fc e ^ is a discrete sequence of parameter dependent eigenvalues corresponding to the 
Dirac eigenvalue problem J^u^x) = X^Uh^x). 

Proof . See [13 E]. ■ 

We will as a motivating example consider perturbations which are locally periodic and of the 
form Vh(x) = V\{x)V2{hx). The entries of V\ are assumed to be (— l)-homogeneous. The entries 
of ^2(2/) are assumed to be periodic with respect to a regular lattice in M 3 . This can also be 
phrased that they are defined on the unit torus T 3 . 

The evolution equation associated with the Dirac operator J^ reads 



(7) > •■" 



u h (-,0)=u£. 



By the Stone theorem, since J%h is self-adjoint on X, there exists a unique solution u^ to (|7|) 
given by 

(8) u h (;t)=U h (t)u° h ,\/u° h eX, 

where U^,(t) = exp(— (i/H)Jtffit) is the strongly continuous unitary operator generated by the 
infinitesimal operator —{i/K)Mh on y, see e.g. [B] or [T4"j . 

In the sequel we will consider a shifted family of Dirac operators denoted by J^ and defined 
as Jffh = H + Vh, where H = H + mc 2 I. Also without loss of generality we will in the sequel 
put h = c = m = 1. By Theorem 1, for any h & N, we then get 

a(J^ h ) = (-00, 0] U {XfokeN U [2, to). 

2.2. G-convergence. For more detailed information on G-convergence we refer to e.g. [H [TT] 
for the application to elliptic and parabolic partial differential operators, and to the monograph 
[3] for the application to general self-adjoint operators. Here we recall some basic facts about 
G-convergence for self-adjoint operators in V. 
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In the present work we frequently write Ah converges to A when we mean that the sequence 
{Ah} converges to A. Let A > 0, by CP^(y ) we denote the class of self-adjoint operators iona 
closed linear subspace "V = T)(A) of V such that (Au,u) > A||w||2 Vit € D(A). 

Definition 1. Let A > 0, and let {Ah} C f\(^) then we say that Ah G-converges to A E J'a(V), 
denoted Ah — : — > A in V if AT PhU — ^-> j4 _1 Pu in ^ Vit G V, where s and w refer to strong 

G,w n w 

and weak topologies respectively, and Ph and P are the orthogonal projections onto % := ~D(Ah) 
and V := T)(A) respectively. Also we say {Ah} C ^(V) converges to {A} C CPo(^) in i/ie strong 
resolvent sense (SRS) if (XI + A fc ) -^-> (A/ + A) in y VA > 0. 

The following result provides a useful criterion for G-convergence of self-adjoint operators. 
See (3| for a proof. 

Lemma 1. Given A > 0, {Ah} C ^(V) and the orthogonal projection Ph onto %. Suppose that 
for every u € V, A^ P/jU converges strongly (resp. weakly) in }$, then there exists an operator 
A G IPaQO suc/i i/tai A — ^-> A ("resp. A G ' w > A) in CP A (y). 

Prom now on we will just use the word "converge" instead of saying "strongly converge", 
hence Ah > A instead of Ah — — > A. 

2.3. Some Basic Results in Spectral Theory. For more details see [21 [U [15]. Given a 

Hilbert space X, let (W,srf) be a measurable space for ^ C C and s$ being a a— algebra on 

fy . Let P = F(X) be the set of orthogonal projections on X, then £7 : &/ — > P is called a 

spectral measure if it satisfies the following 

(i) .£7(0) = (This condition is superfluous given the next properties). 

(ii) Completeness; E(fy) = I. 

(in) Countable additivity; if {A n } C stf is a finite or a countable set of disjoint elements and 

A= U n A„, then E(a) = ^ E(A n ). 

n 

If £7 is spectral measure then E(A± n A2) = £7(Ai)£7(A2) = £7(A2)£7(Ai), also £7 is modular, 

i.e. £7(Ai U A2) + £7(Ai n A2) = £7(Ai) + £7(A2). For an increasing sequence of sets A n , 

lim £7(A n ) = £7(U n A„), while if A n is a decreasing sequence then lim E(A n ) = E(n n A n ). 

n—toc ' n—too 

Because of the idempotence property of the spectral measure we have ||£7 n ?x|| 2 = (£7 n u, £7„u) = 
{E^u, u) = (E n u, u) — > (Eu, u) = ||£7u|| 2 , which means that the weak convergence and the strong 
convergence of a sequence of spectral measures £7 n are equivalent. 
Let E u ^ u (a) be the finite scalar measure on srf generated by £7, 

E U , U (A) = (E(A)u,u) = \\E(A)u\\ 2 

and E U)V (a) be the complex measure 

£7 U ,„(A) = {E(A)u, v),Vu,v€X. 
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By the above notations E„ v ( A) < ||iiYA)ii||||.E(A)u|| < ||«|||M|. 

Let *% = R. The spectral measure on the real line corresponding to an operator S is denoted 
by £/ (A) (where the superscript S indicates that the spectral measure E corresponds to a 
specific operator S) 

E s (\) = E s (a) , where A= (-00, A) , for A € M . 

Clearly E (A) ismonotonic (nondecreasing) , i.e. E (Ai) < E (A2) for Ai < A2. Also lim E (A) 

A— >— 00 

and lim E (A) = I. E S (X) is self-adjoint, idempotent, positive, bounded, right continuous 

A— >oo 

operator ( lim E (A + t) = E (A)), and discontinuous at each eigenvalue of the spectrum. If A 

is an eigenvalue, then we define p(X) = E S (X) — E S (X — 0) to be the point projection onto the 
eigenspace of A. For A being in the continuous spectrum p(X) = 0. 
Now we state the spectral theorem for self-adjoint operators. 

Theorem 2. For a self-adjoint operator S defined on a Hilbert space X there exists a unique 
spectral measure E s on X such that 

(i) S= f XdE s (X). 

Ja(S) 

(») E{A) =0 if Ano-(S) =0. 
(Hi) I/ACK is open and ADa(S) ^ 0, then E(a) / 0. 

Proof . See e.g. [2J . ■ 

3. The main results 

Consider the family {Jifh}h^m of Dirac operators with domain D(J%) = X. We will state and 
prove some useful theorems for operators of the class ^xi^) for A > 0, where X and V are the 
Hilbert spaces defined above. The theorems are valid for general Hilbert spaces. 

The following theorem gives a bound for the inverse of operators of the class 7\(^) for A > 0. 

Theorem 3. Let A be a positive and self-adjoint operator on y and put B = A + XI. Then for 
A > 

(i) B is injective. Moreover, for every v £ H(B), (B~ 1 v,v) > A||-B -1 f||u and \\B~ 1 v\\y < 

A" 1 HI*. 

(**) R(#) = y. 

Proof. See Propositions 12.1 and 12.3 in |3j. 

The connection between the eigenvalue problems of the operator and its G-limit of the class 
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"PxCd) for A > is addressed in the next two theorems. Here we prove the critical case when 
A = 0, where for A > the proof is analogous and even simpler. 

Theorem 4. Given a family of operators {Ah} of the class CPo(y) G-converging to A G CPo(^) 
in the strong resolvent sense. Let Uh be the solution of A^u^ = fh, where {fh} is converging to 
f in]j. If {uh} converges to u in^, then u solves the G-limit problem Au = /. 

Proof . Since Ah G-converges to A in the strong resolvent sense 

(9) B- x P h v -> B- l Pv , W G V , 

where Bh and B are Ah + XI and A + XI respectively. Note that by Theorem 4, T)(BT ) = 
R(Bh) = V, so the projections P re and P are unnecessary. 

Consider A^Uh = fh which is equivalent to BhUh = fh + Att ra , by the definition of Bh we have 
Uh = B^ (fh + Xuh)- Define J?h = fh + Aii ra which is clearly convergent to J* = f + Xu in y by 
the assumptions. Therefore -B^ ^^ — ► S _1 # , this is because 

< ll^lb II A - ,/lla + KV - B^fh 
— >o. 

The convergence to zero follows with help of ^ and the boundedness of the inverse operator 
B h l . Thus, for all v G y 

(«,«) = lim (« h ,u) = lim {B^^v) = {B' 1 J » . 

h^tca re— foo 

Hence (u — B^ 1 ^ ,v) = for every uG^, which implies -Bit = ^, therefore ^4u = /. ■ 

Theorem 5. Ze£ {^4^} 6e a sequence in ^(V) which G-converges to A G !Po(y) in i/ie strong 
resolvent sense, and let {^h, u h} be the solution of the eigenvalue problem AhUh = fJ-h u h- If 
{/J,h, Uh} ~> {m? n l ™Kxy, i/ien i/ie limit couple {fi, u} is the solution of the eigenvalue problem 
Au = jiu. 

Proof . The proof is straight forward by assuming /„ = fih u h (which converges to \xu in y) in 
the previous theorem. ■ 

The convergence properties of self-adjoint operators has quite different implications on the 
asymptotic behavior of the spectrum, in particular on the asymptotic behavior of the eigenval- 
ues, depending on the type of convergence. For a sequence {Ah} of operators which converges 
uniformly to a limit operator A nice results can be drawn for the spectrum. Exactly speaking 
{a (Ah)} converges to cr(A) including the isolated eigenvalues. The same conclusion holds if the 
uniform convergence is replaced by the uniform resolvent convergence, see e.g. [6] . In the case of 
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strong convergence (the same for strong resolvent convergence), if the sequence {A^} is strongly 
convergent to A, then every A G o~(A) is the limit of a sequence {A/j} where A^ G a(Ah), but 
not the limit of every such sequence {A^} lies in the spectrum of A, (see the below example 
taken from [E]). For weakly convergent sequences of operators no spectral implications can be 
extracted. 

Example. Let A^ be an operator in L 2 (R) defined by 

d 2 
A i>h = -—j + V it h(x) , for h G N and i = 1, 2 , 

where 

/ "!» if/i<^</i + l, , T/ , . J -1, ifx>/i, 

riM^) = s an d ^^(x) = < 

I , Otherwise , I , Otherwise . 

The operator A^ converges to A = —4-^ in the strong resolvent sense as h — > oo for both % = 
1, 2. One can compute the spectrum for the three operators and obtain a(Ai < f l ) = [0, oo) U {nh} 
for /ih being a simple eigenvalue in [—1, 0] and o-(A 2! h) = [—1, oo), whereas for the unperturbed 
limit operator A the spectrum consists of just the continuous spectrum, i.e. cr(A) = [0, oo). 

Since the uniform convergence is not always the case for operators, the theorem below provides 
some criteria by which the G-convergence of an operator in the set 3^(9) (and hence the G- 
convergence in the strong resolvent sense of operators of the class CPo(^ )) implies the convergence 
of the corresponding eigenvalues. 

Theorem 6. Let {A^} be a family of operators in l ?x(^), A > ; with domain X. If A^ G- 
converges to A G ^xi^), then %h '■= A^ converges in the norm of £(V) (tB(V) is the set of 
bounded linear operators on ^) to % := A~ l . Moreover the k eigenvalue n\ of Ah converges 
to the k eigenvalue fj, of A and the associated k eigenvector u^ converges to u weakly in X, 
\/k GN. 

Proof . By the definition of supremum norm 

(10) \\0Cfi — X||$(y) = sup \\XhV — 0Cv\\y = sup \\%hV — "Kv\\y . 

IMIy =1 Nly<i 

Also, by the definition of supremum norm there exists a sequence {vh} C ^ with ||fh||y < 1 such 
that 

(11) \\X h - 3C|| S (« < \\X h v h - Xv h \\y + I ■ 

It is well-known that %h and % are compact self-adjoint operators on y. Both are bounded 
operators, by Theorem 3, with compact range X of y. 
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Consider now the right hand side of (jlip . We write this as 

\\X h V h - %V h \\y + -< \\X h V h - X h v\\y + \\X h V - Xv\\y + \\Xv h - Xv\\y + -. 

The first and the third terms converge to zero by the compactness of Xh and X on y and the 
second term converges to zero by the G-convergence of Ah to A. Consequently 

(12) ||3Cfc - 3Cb(v) ->• 0. 
Consider the eigenvalue problems associated to A^ and A^ 1 

(13) A-'vt = X k h v k h ,ken 
and 

(14) A~ 1 v k = X k v k , ken. 

Since A^ and A^ 1 are compact and self-adjoint operators it is well-known that there exists 
infinite sequences of eigenvalues A^ > X\ > ■ ■ ■ and A 1 > A 2 > • • • accumulating at the origin, 
respectively. Define pJi := (A^) _1 and /i fc := (A fe ) _1 for all k £ N. Consider now the spectral 
problems associated to Ah and A 

(15) A h u\ = fx k h u k h , k € N 
and 

(16) Au k = n k u k , ken. 

There exists infinite sequences of eigenvalues < n\ < n\ < ■ ■ ■ and < fi 1 < [i 2 < ■ ■ ■ 
respectively. By the compactness of Xh and X the sets {\ k l }'^ =1 and {X k }^Li are bounded in M, 
thus the proof is complete by virtue of the following lemma. ■ 

Lemma 2. Let X, y, Xh, X, X\ and A be as in Theorem 6, and let Ah € J'a(V), A > 0. There 
is a sequence {rfy converging to zero with < r^ < A such that 



\k 

(17) |A^ - A fc | < c— E sup \\X h u-Xuh, 

A ~ r h u£^(\ h ,X) 

IMIy=l 



where c is a constant independent of h, and ^Y(X k ,X) = {u € D(X) ; Xu = X k u} is the 
eigenspace of X corresponding to X k . 

Proof. See Theorem 1.4 and Lemma 1.6 in [8] Chapter 3. ■ 
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We can now complete the proof of Theorem 6. By the G-convergence of A^ to A we obtain, 
by using Lemma 2 and f)12[) . convergence of the eigenvalues and eigenvectors, i.e. [x^ — >■ /i fc and 
u k _£, u k wea kiy [ n x as h — > oo. ■ 

Let us now return to the shifted and perturbed Dirac operator .3^h- We will throughout this 
section assume the hypotheses of Theorem 1. We further assume that the 4x4 matrix- valued 
function Vh is of the form Vh(x) = Vi(x)V2(hx) where V\ is (-l)-homogeneous and where the 
entries of V<i{y) are 1-periodic in y, i.e. 

v* j (v + k) = vZ\v), fcez 3 . 

We also assume that the entries of V2 belong to L°°(M 3 ). It is then well-known that 



(18) V?{hx) -)■ M(V?) = / Vi j (y)dy, 



in L°°(M 3 ) weakly*, where T 3 is the unit torus in M 3 . It easily also follows from this mean- value 
property that 

V h -> ViM(V 2 ), 

in L P (R 3 ) weakly for p > 3, cf the hypotheses in Theorem 1. 



We are now interested in the asymptotic behavior of the operator and the spectrum of the 
perturbed Dirac operator J^%. We recall the spectral problem for J%, i.e. 

J^ h u h (x) = X h u h {x) 

where there exists a discrete set of eigenvalues {A^}, k = 1, 2, ... and a corresponding set 
of mutually orthogonal eigenfunctions {u^}. We know, by Theorem 1, that the eigenvalues 
(or point spectrum) a p {Mh) C (0,2). We also know that J^ has a continuous spectrum 
a c ( r Jtf?h) = (— 00, 0] U [2, 00). This means that the Dirac operator is neither a positive or negative 
(semi-definite) operator and thus the G-convergence method introduced in the previous sec- 
tion for positive self-adjoint operators is not directly applicable. In order to use G-convergence 
methods for the asymptotic analysis of J#h we therefore use spectral projection and study the 
corresponding asymptotic behavior of projections J£/j which are positive so that G-convergence 
methods apply. 

Let &/ be a fixed c-algebra of subsets of M, and let (M, &/) be a measurable space. Consider the 
spectral measures E and E h of the families of Dirac operators Mh and H respectively, each 
one of these measures maps si onto P , where P is the set of orthogonal projections on X. By 
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the spectral theorem 

(19) Mk= f XdE^ h (X). 
By the spectral theorem we can also write 

(20) f XdE*(X) + V h , 

Ja(U) 

since Vh is a multiplication operator. 

We recall that D(J%) = X, let now 

,yV^ = {u h € X; je h u h = X k h u h }, 

i.e. the eigenspace of J&h corresponding to the eigenvalue X\. Further define the sum of mutual 
orthogonal eigenspaces 

where X p h is a closed subspace of ^ invariant with respect to J^,. 



It is clear that for Uh S X p h we have 

(J$? h u h ,u h ) = X k \u h \ 2 > 0, A; = 1,2,.... 
Let us now consider the restriction r J%? p of J4fh to X p h which can be written as 

3^= Y. A^^(A), 

where the spectral measure E h ' p is the point measure, i.e. the orthogonal projection onto 
ker(J4?h — XI). With this construction J4? p is a positive and self-adjoint operator on X p h with 
compact inverse {J^ p )~ l . By Lemma 1, see also Proposition 13.4 in [3], we conclude that there 
exists a positive and self-adjoint operator ^C v such that, up to a subsequence, M? p G-converges 
to J& p , where JT P has domain D(^ p ) = X p where 

X p = e fceN ^ /fc 

is a closed subspace of y and where 

J^ k = {ueX;J^ p u = X k u}. 

Moreover, by Theorem 6, the sequence of k th eigenvalues {A^} associated to the sequence {jifff } 
converges to the k eigenvalue of X k of J^f p and the corresponding sequence {«<!} converges to 
u k weakly in X. The limit shifted Dirac operator restricted to X p is explicitly given by 

Jg> p = (k + V 1 M(V 2 ))\ X p- 
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This follows by standard arguments in homogenization theory, see e.g. [T]. 

We continue now to study the asymptotic analysis of the projection to the closed subspace 
of y corresponding to the positive part [2, +00) of the continuous spectrum of ,!%h- 

We denote by X c h the orthogonal complement in X to the eigenspace X p h . Thus, X c h is the 
closed subspace invariant with respect to M^ corresponding to the absolutely continuous spec- 
trum a c (Mh) = (—oo,0] U [2,oo). We now define the two mutually orthogonal subspaces X^ 
and X c £~ with 






where XJ corresponds to the positive part [2, +cx>) and X^ corresponds to the negative part 

c " 



(—oo,0], respectively. Next we define the restriction Jtff' of J%h to X 1 ? which can be written 



as 



^ c ' + = f XdE^' c > + (X), 

where the spectral measure E h ' c,+ (X) is the continuous spectral measure corresponding to 
Jtf£' . By construction J^' is a positive and self-adjoint operator on X?' . Therefore by 
Proposition 13.4 in [3], there exists a sequence {J^' } which G-converges to a positive and 
self-adjoint operator J4? c ' + € X c ' + . Moreover, since A is not an eigenvalue, the corresponding 
sequence {E h ' c,+ (X)} of spectral measures converges to the spectral measure E ,C ' + {X) corre- 
sponding to J^ c,+ . 



Let us consider the evolution equation 
(21) 



i§- t u. h (t,x) = J^' + u h (t,x) 



{ Ufc(-,0)=u°. 

By the Stone theorem, there exists a unique solution u^ = u(x, t) to ([2T]) given by 

U h (;t)=U h (t)u° h , Vu° G X£ + , 

where U/ l (t) = exp(— iJ^' t) is the strongly continuous unitary group of transformations gener- 
ated by the infinitesimal operator — iffl^ on y. By the G-convergence of the sequence {J^^ } 
it follows that the associated sequence {1X^' (£)} of unitary groups of transformations converges 
to a unitary group of transformations U c ' + (t) which for every u° G X c,+ defines the solution 
u(-,t) = lt(t)u° to the limit evolution equation 

j ifu ( t,x) = J^ c ' + u(t,x), 
1 u(-,0)=u°. 
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Finally, by considering the operator —fflu 



where M 3 ^' 



is the restriction to X?' , i.e. the 



closed subspace corresponding to the negative part (—00, 0] of the continuous spectrum we can 
repeat all the arguments from the positive part of the continuous spectrum. 
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